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^ ■ Abstract 

Qh| Motivated by the possible large direct CP asymmetry of B^ — ► D + D decay measured 

by Belle collaboration, we investigate double charm B u ^ and -B s decays in the minimal 
\ supersymmetric standard model with R-parity violation. We derive the bounds on rele- 

\q • vant R-parity violating couplings from the current experimental data, which show quite 



consistent measurements among relative collaborations. Using the constrained param- 
eter spaces, we explore R-parity violating effects on other observables in these decays, 
which have not been measured or have not been well measured yet. We find that the 
R-parity violating effects on the mixing-induced CP asymmetries of — > D* + D* 
and B® — > D^* + decays could be very large, nevertheless the R-parity violating ef- 

X 

fects on the direct CP asymmetries could not be large enough to explain the large direct 
CP violation of B% — > D D from Belle. Our results could be used to probe R-parity 
violating effects and will correlate with searches for direct R-parity violating signals in 
future experiments. 
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1 Introduction 



Double charm decays of B u ^ and B s provide us with a rich field to study CP violation and 
final-state interactions as well as to extract information of Cabibbo-Kobayashi-Maskawa (CKM) 
elements. CP asymmetries (CPAs) in these decays play important roles in testing the Standard 
Model (SM) as well as exploring new physics (NP) 

Double charm decays, B® — > D { * )+ D { * ] , B~ — > D { * )0 D { * ] and B® — > D* )+ D* , are 
dominated by color-allowed tree 6 — » ceo? transition, but involve small penguin pollution from 
b — > nuci transition carrying a different weak phase. The latter contributions lead to direct 
CPAs, which are very small (about the order of 10~ 2 ) in the SM. If penguin corrections are 
neglected, the SM predictions for the direct CPAs would be zero. It is interesting to note that 
both BABAR and Belle have measured the direct CPA in B® — > D + D~ decay 



-0.91 ± 0.23 ± 0.06 (Belle [3j), 
-0.07 ± 0.23 ± 0.03 (BABAR [4]), 



C(BIB»^D + D-) = 
respectively. One would find the difference between the two measurements is 



(1) 



AC = 0.84 ±0.32, 



(2) 



i.e, the difference is as large as 2.7a. So far, such a large direct CPA has not been observed 



in the other measurements of B® 
Q, 5 , 1, E , I, Q, lo], 11 ] , which i 



D M+ D M -, B- 



D ( * )0 D M - and B° 



£)(*)+£)(*) 



involve the same quark level weak decays. If the large CP violation 
in Bj — > D + D~ from Belle is true, it would establish the presence of NP. At present one cannot 
conclude the presence of NP in those decays. Equivalently one also cannot take the CPAs are 
in agreement with the SM expectations. Recently the large direct CPA in _B9 — s- D + D~ has 



been investigated with possible NP scenarios, such as unparticle interaction 12| and the NP 
effects in electroweak penguin sector [k], [ljj] , and so on. 

In this paper, we would like to investigate B u ^ and B„ double charm decays systematically in 



the minimal supersymmetric standard model (MSSM) 15 



16[ with R-parity violation 



In the literature, the possible appearance of the R-parity violating (RPV) couplings 



17 



17 



which violate the lepton and/or 



searching for supersymmetry 19 



Daryon number conservations, has gained full attention in 



20| | . The effects of supersymmetry with R- 



jarity vio 



in B meson decays have been extensively investigated, for instance in Refs 
In our work, twenty-four double charm decays B^ 



21 



22 



23 



ation 



m. 



D M+ D M -, B- 



D ( * )0 DZ and B° s 



25] 



are studied in the RPV MSSM. For simplicity we employ naive factorization 
for the hadronic dynamics, which is expected to be reliable for the color-allowed amplitudes, 
which are dominant contributions in those double charm decays. 

The color-allowed tree level dominated decays of b — > ccd, i.e. B® — > D* + D* , B~ — > 



D M °D M - and B? 



10 



11 



26|]. 



involve the same set of RPV coupling constants. For these 
processes, besides the CPA in B® — > D + D~, a few other observables in Bj. d -> D^DM have 
been already measured by BABAR and Belle collaborations {3], 4, 0, O, Q, 0, 
To derive constraints on the relevant RPV couplings, we will choose a set of data from the 
aforementioned measurements which have quite high consistency between the measurements of 
BABAR and Belle. Then, using the constrained RPV coupling parameter spaces, we predict 

B - _> D W0 D W ~ and B°„ -> 



the RPV effects on the other observables in B® 



D { * )+ D (t) decays, whose measurements from BABAR and Belle are not compatible within 2a 
error range, and/or which have not been measured yet. One of our goals is to see how large the 
direct CPA of B d — > D + D~ can be within the constrained parameter spaces. We find that the 
lower limit of C(B%, B® — > D + D~) could be just slightly decreased by the RPV couplings, and 
the RPV effects on this quantity are not large enough to explain the large direct CP violation 
from Belle, although the mixing -induced CPAs of B° d -> D D are very sensitive to the 



, B- 



D <t)0 D { ; ] - and 5° 



D* + D* are governed by the b 



RPV couplings. 

Decays B° d -> D w+ D™~ 
ccs transition at the quark level, and also involve the same set of RPV coupling constants. They 
have similar properties to B® — > D ( * >+ D { * } , B~ — > D*°D* and B® — > D { * )+ D { * ] decays, 



nevertheless the penguin effects are less Cabibbo-suppressed. For these 



decays, most branching 



ratios and one longitudinal polarization have been measured [27], [28|, [29|, [30|, |3ll, |32|, [33j, [34|. We 
will take the same strategy as the one for b — > ccd decays to constrain relevant RPV couplings 
and estimate RPV effects in these decays. We find that RPV couplings could significantly affect 
the CPAs of these decays, and could flip their signs. 

Our paper is organized as follows: In Sec. 2, we briefly introduce the theoretical framework 
for the double charm B u d and B s decays in the RPV MSSM, and we tabulate all the theoretical 
input parameters. In Sec. 3, we deal with the numerical results and our discussions. At first, 
we give the SM predictions with full uncertainties of the input parameters. Then, we derive 
the constrained parameter spaces which satisfy all the experimental data with high consistency 
between different collaborations. Finally, we predict the RPV effects on other quantities, which 



3 



have not been measured or have not been well measured yet. Section 4 contains our summary. 



2 Theoretical Framework 



2.1 Decay amplitudes in the SM 

n the SM, the low energy effective Hamiltonian for AB = 1 transition at a scale \x is given by 
35] 



n !s = -t= E K\ CxQl + c 2 q\ + £ c< + cv 7 q 77 + c 8 ,q 8 , \+h.c, (3) 

V ^ p=u,c I i=3 ) 

here A p = V p bV* for 6 — > g transition (p G {u,c},g G {d,s}). The detailed definition of the 



effective Hamiltonian can be found in 



It is empirically observed that naive factorization 25[ still works reasonably well in the 
color-allowed double charm B u ^ and B s decay processes. We will describe the B — > D* D^ 1 
decay amplitudes within the naive factorization approximation in this paper. Under the naive 
factorization approximation, the factorized matrix elements are given by 



.4 



[BD i * ) ,D ( * ) ] 



D ( ; } \q^(l - 75 )c|0) (D M |c Tm (1 - l5 )b\B) . 



(4) 



Decay constants and form factors |36|, 37] are usually defined as 



(D;(p D ,)\qYc\0) 

H 9 

(D(p D )\c^b\B(p B )) 
(D*(p D ^e*)\c^b\B(p B )) 
(D*(p D *,e*)\c W b\B(p B )) 



2V(q 2 ^ 



m-r, — m 
a d 



(5) 
(6) 



m B + m D , 



(8) 



e*{m B + m Dt )A 1 (q 2 ) - (p B +p D ,) li (e* ■ Pi 



A 2 (q 2 



m B + m D *_ 



2m, 



■[A 3 (q 2 ) - A (q% 



(9) 



with q = P b — P D (*)- Then we can express A 



BD (*) D (*)i in terms of decay constants and form 



4 



factors as follows 



j4 [BO (, »,flW] 



ifD q (m 2 B -m 2 D )F Q (m 2 Di 
2fD* m B\Pc\Fi(m 2 D .) 

i q 

-2f Dq m B \p c \Ao(m 2 Dn ) 
-if D *m 



(DD* q ), 

(D*D q ), 



1 Q 



){m B +m D *)A 1 {m 2 D *) 
i i 

-K* ■ Pd;)K ; ■ Pd*)^^^ 

2V(m 2 * ) 
i ■ -^/^ *u <y. ft ®q 

-rt£fiva{3£ D *£ D *P D *P D * ms+m 

q q D D 



(10) 



[D*D*). 



Decays B — > D * D* may occur through both tree level and loop induced (penguin) quark 
diagrams, and the SM decay amplitudes within the naive factorization are given as 



M SM (B - D W D^) = ^ [\ c a\ + J^X, [a\ + a% + tfog + «!)]) (11) 



where the coefficients <£ 



Pf with the upper (lower) sign applied when i is odd 



(even), and P[ account for penguin contractions. The factorization parameter £ in Eq. (fTTl) 
arises from the transformation of (V — A)(V + A) currents into (V — A)(V — A) ones for the 
penguin operators Qs, ■ • ■ , Q$, and it depends on properties of the final-state mesons 



2m- r 



(m c +m q )(m b -m c ) 



2m lq 

(m c +m q )(m b +m c ) 





(DD g ), 
{DD* q ), 
(D*D q ), 
(D*D*). 



(12) 



For the penguin contractions, we will consider not only QCD and electroweak penguin opera- 
tor contributions but also contributions from the electromagnetic and chromomagnetic dipole 
operators. Pf are given as follows 



Pi 



pp 



pP 
^8 



pP 



G D w(mp) 



0, 

9vr I L 9 

— Ic \— 

9vr I 1 [ 9 



G (*)(m p ) 



10 
~9~ 
10 







-2F 2 C|f] 




- G (*)(m p ) 


- 3F 2 C% f 


- G (*)(m p ) 





(13) 
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where the penguin loop-integral function G (,){m p ) is given by 



D 



duG(m„, fc)<3> (*) (u), 



G(m p , k) 



—4 / dxxil — x)hi 
Jo 



m 2 — k 2 x(l — x) 
ml 



— te 



(14) 
(15) 



with the penguin momentum transfer k 2 = m 2 c + u(m 2 — m 2 — m? M ) + u 2 m 2 M2 , where u = 
1 — u. In the function G n (*)(m p ), we have used a I)W meson-emitting distribution amplitude 



$ D (*)(n) = 6u(l —«)[! + a n (*)(l — 2u)], in stead of keeping k 2 as a free parameter as usual. 



The constants F\ and F2 in Eq. ( fl~3l) are defined by 



ml-um 2 Di] -2ml+m b mc 



m b —m c 



(DD. 



J du<S> D .{v)*$ [um b + 



2umj 



m b -m c 



-e* 2 ■ pi — um c 



m^—um 2 ^ —2m 2 —rn b in c 

1 v ™J m b +m c k 2 
2um T 



(D*Dg), 

(D*D*), 



where e\ L • p\ 



[mt 



fodu$ Dq (u 

/ du$ D *(u)?§t (um b + ^t/ 2 ■ Pi + um c 

Jo 1 du<$> Dq (u)^[u(m b - m c ) + m c ] (DD q ), 

(DD*), 

Jo 1 du$ Dq {u)^[u{m b + m c ) - m c ) (D*D q ), 

(D*D*), 

m 2 M * - m 2 c ) / (2m M *) and e* 2T ■ p 1 = for B —> D*D* decays. 



(16) 



(17) 



2.2 Decay amplitudes of the RPV contributions 



In the RPV MSSM, in terms of the RPV superpotential [13], we can obtain the relative RPV 
effective Hamiltonian for B — > D ( * ] ' D { * ] decays as following 



n/KPV 



E 



A ikn A jln rn -4 : /l3 

2m 2 , 



2 V 



■{dkYPRU j ) 1 {ui'y IJ ,P R di) 1 + (dk^PRU^siua^PRdijs 



\i y* 

* nip ' * -iV 



+ E li2 ^^- 8//30 (4PL^)l(«„P^)l + 



(18) 



where Pl 



^2^,1] = and /3 = 11 — f^/- The subscripts 1 and 8 of the 



currents represent the currents in the color singlet and octet, respectively. The coefficients 
yj-4/Po anc i ^-8/ A) are dug to the running from sfermion mass scale mj (assumed as 100 GeV) 
down to m b scale. Since it is usually assumed in phenomenology for numerical display that 
only one sfermion contributes at one time, we neglect the mixing between the operators when 
we use the renormalization group equation (RGE) to run 7"^f v down to the low scale. 
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The decay amplitudes of RPV contributions to B — > D ( * ] D l * ] are given by 



A" 



•1 + ^) (D? |?7m(1 + 7s)c|0) (£> W |cTm(1 + ^)b\B 



2A' (^f - 75)c|0) (D w |c(l + 75 )6|5 



^(*) I 



-A"(-l + £)+£A' 



.4 



[BD,D 



A " I" 1 + ife) 



,4 



[BD,D*] 



A" ( -1 + ifc ) - £A' 



A " + it) 



.4 



[BD*,D q 



[BD*,D*] 



(DO;), 

(D*D q ), 

(D*D*), 



(19) 



(20) 



where A" = r)~ 4 ^ A ^| 12 ( ^f^ ) and A' = r]~ 8 ^ Ei ^f^ 21 ( A f# 22 ) for g = d (g = s). 
^'[_bd*,d*] i s defined by 



A '[BD',D' q ] = % fo* q m D q 

-( £ *d* -Pd*)(£*d* -Pd 



(e* D * ■ e* D *)(m B + m D *)Ai(m^*) 

9 9 



2A 2 (m 2 D * 



2V(m 2 «) 

1 ; ^nuapt n* £ D*PD*Pn* I 



.(21) 



2.3 Observables to be investigated 

We can get the total decay amplitudes in the RPV MSSM as 

M{B - D ( * ] D ( * ] ) = M SM (B - D ( * ] D ( * ] ) + .M RPV (5 - D W D^). 
The branching ratio *B reads as 



(22) 



-6(5 -> D M D M ] 



]irm 2 B 



(*) 



(23) 



where tb is the i? lifetime, \p c \ is the center of mass momentum in the center of mass frame of 
B meson. In B — > D*D* decays, the two vector mesons have the same helicity, therefore three 
different polarization states, one longitudinal and two transverse, are possible. We define the 
corresponding amplitudes as -Mq,± m the helicity basis and A4l,\\,± m the transversity basis, 



which are related by .Ml = -Mo and Ai 



M+±M- 
V2 



Then we have 



m(b -»• d*d: 



\M \ 2 + \M + \ 2 + \M^\ 2 = \M L \ 2 + \M {l \ 2 + \M ± \ 2 . 



(24) 



The longitudinal polarization fraction fi and transverse polarization fraction f± are defined by 



f Lt± (B - d*d; 





M L ,i 


2 


\M l 


2 + 


\M l{ \ 


2 + \M ± \ 2 



(25) 
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In charged B meson decays, where mixing effects are absent, the only possible source of 
CPAs is 



A 



fc,dir 
CP 



M k (B- - f)/M k (B+ - 


-/) 


2 

- 1 


M k (B- - J)/M k (B+ - 


-/) 


2 

+ 1 



(26) 



and k 



for B~ -> D*D* decays and k — L for 5" -> DD,, DD* D*D q decays. Then 



for B u — ► D*D* decays, we have 



4fl-M|| 


2 + A L c f\M L 


2 




M\\ 


2 +\M L 


2 



(27) 



For CPAs of neutral B q meson decays, there is an additional complication due to B q — B q 
mixing. There are four cases that one encounters for neutral B q decays, as discussed in Refs. 



38 



39 



40 



4l|. 



i) B® — > f, B® — > f, where / or / is not a common final state of B q and B®, for example 



q 

D+D 



(ii) B® —►(/ = /) <— B® with f CP = ±/, involving final states which are CP eigenstates, 
i.e., decays such as B® — > D + D~ , B® — > D+Dj . 

[in) B® —> (f = f) <— B® with / CP ^ ±/, involving final states which are not CP eigenstates. 
They include decays such as B q — > (VV)°, as the VV states are not CP eigenstates. 



(iv) B° q 



— > (/&/) <— with / ^ f, i.e., both / and / are common final states of 
B° q and B°, but they are not CP eigenstates. Decays B Q d (B° d ) -> D*-D + ,D~D* + and 
B°(-B°) -> D*~D+,D~D* + belong to this case. 

CPAs of neutral B decays in case (i) are similar to CPAs of the charged B decays, and 
there are only direct CPAs A^p since no mixing is involved for these decays. For cases (ii) 
and (iii), their CPAs would involve B q — B q mixing. The time-dependent asymmetries can be 
conveniently expressed as 



A){t) = S k f sin(Ami) - C) cos(Amt), 



S) 



2Im(A jS 
1 + lA, 



|2 ' 



nk 



|2 ' 



(28) 
(29) 



where At = a !1 fe lno f ■ In addition, <St and Ct can be obtained from the similar relation given 
in Eq. ([57]). 



Case (iv) also involves mixing but requires additional formulae. Here one studies the four 



time-dependent decay widths for B°(t) -> /, B°(t) -> /, B°(t) -> / and 5°(t) -> / 



411 ] . These time-dependent widths can be expressed by four basic matrix elements 



138 

3 



39 



4d, 



s f = 


2Im(A' fe ) 


nk 


l - 


|A' fc 


2 


1 + KI 2 ' 


l + 


|A' fc 


2 • 


qh _ 


2Im(A») 


nk 


l - 


|A2 


2 


1 + 


l + 


|A2 


2 ' 



g = (f\H eff \B° q ), h=(f\H eff \B° q ), 

g = (f\Heff\B° q ), h=(f\H eff \B° g ), (30) 

which determine the decay matrix elements of B q —+ f,f and of B q — > f, f at t — 0. We will 
study the following quantities 

(31) 
(32) 

with X' k = (q/p)(h/g) and = (q/p)(g/h). The signatures of CP violation are T(Bg(t) —> 
/) ^ r(5°(t) -> /) and T(B° q (t) -»• /) ^ r(B°(t) -> /), which means that C f ^ -C/ and/or 
5/ ^ -5/. 

2.4 Input parameters 

Theoretical input parameters are collected in Table [U In our numerical results, we will use the 
input parameters which are varied randomly within la range. 
We have several remarks on the input parameters: 

• CKM matrix elements : The weak phase 7 is well constrained in the SM, however, with 
the presence of R-parity violation, this constraint may be relaxed. We will not take 7 
within the SM range, but vary it randomly in the range of to n to obtain conservative 
limits on RPV couplings. 

• Decay constants: The decay constants of D* mesons have not been directly measured in 
experiments so far. In the heavy-quark limit (m c — > 00), spin symmetry predicts that 
Id* = /d„, and most theoretical predictions indicate that symmetry-breaking corrections 



45 



46|. Hence, we take f D . = (1.1 - 1.2) f Dq 



enhance the ratio fo* q l Jo q by 10% — 20% 
as our input values. 

Distribution amplitudes: The distribution amplitudes of D q *' mesons are less constrained, 
and we use the shape parameter a D {*) = 0.7 ± 0.2 and a n (») = 0.3 ± 0.2. 



Table 1: Summary of theoretical input parameters and ±lcr error ranges for sensitive parameters 
used in our numerical calculations. 



m Bu = 5.279 GeV, m Bd = 5.280 GeV, m Bs = 5.366 GeV, M D o = 1.865 GeV, 
M D+ = 1.870 GeV, M n+ = 1.969 GeV, M D *o = 2.007 GeV, M D ,+ = 2.010 GeV, 
M D *+ = 2.107 GeV, m b (m b ) = (4.20 ± 0.07) GeV, m c (m c ) = (1.25 ± 0.09) GeV, 
fn s (2 GeV) = (0.095 ± 0.025) GeV, m u (2 GeV) = (0.0015 ~ 0.0030) GeV, 
m d (2 GeV) = (0.003 ~ 0.007) GeV, 

r Bu = (1.638 ± 0.011) ps, r Bd = (1.530 ± 0.009) ps, r Bs = (1.425lg;g£) ps. [42] 
\V,J = 0.97430 ±0.00019, \V,J = IKJ = 0.00344+ - 00022 



0.00082; I v ub\ ~ u - uuo ^-0.00017) 

cd\ = 0.22508io;ooo82) lKs| = 0.97350lo;oool2; |K&| = 0.04045^o;ooo7g, 
«| = 0.00841+Hg, |^| = 0.03972t°;K, \V tb \ = 0.999176l~, 



90.7^)°, /3=(21.7± o: °)°, 7 =(67.6ll;g)°. \M 

f D = (0.201 ± 0.003 ± 0.017) GeV, f Ds = (0.249 ± 0.003 ± 0.016) GeV. [44J 



Form factors : For the form factors involving B — > transitions, we take expressions 
which include perturbative QCD corrections induced by hard gluon vertex corrections of 
b — > c transitions and power corrections in orders of l/m b c [37], As for Isgur-Wise 



function we use the fit result = 1 - 1.22(w - 1) + 0.85(cu - l) 2 from Ref. [481]. 
Wilson coefficients : We obtain Wilson coefficients in terms of the expressions in 35]. 



RPV couplings: When we study the RPV effects, we consider only one RPV coupling 
product contributes at one time, neglecting the interferences between different RPV cou- 
pling products, but keeping their interferences with the SM amplitude. We assume the 
masses of sfermion are 100 GeV. For other values of the sfermion masses, the bounds on the 
couplings in this paper can be easily obtained by scaling them by factor f 2 = ( loo^Gev ) 2- 



3 Numerical results and discussions 

In this section we summarize our numerical results and analysis in the exclusive color-allowed 
b — > ccq decays. First, we will show our estimates in the SM with full theoretical uncertainties of 
sensitive parameters. Then, we will investigate the RPV effects in the decays. We will constrain 
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relevant RPV couplings only from quite highly consistent experimental data and show the RPV 
MSSM predictions for the other observables, which have not been measured yet or have less 
consistency among different collaborations. 

3.1 Exclusive color-allowed b —> ccd decays 

Decays B d — > D ( * )+ D { * ] , B~ — > D ( * )0 D ( * } and B® — > D ( s * )+ D M are dominated by the color- 
allowed b — > ccd tree diagram, but involve small penguin pollution from the b — * uud transition 
carrying a different weak phase. These decays involve the same set of RPV coupling constants 
A232A212 an d A / i23 A / j 21 at tree level due to squark and slepton exchanges, respectively. For 
B° d -> D w+ D w ~ and B~ -> D w ' processes, a few observables have been measured 
by BABAR and Belle collaborations. The latest experimental data and their weight averages 
are summarized in Table EJ We can see almost all physical quantities have been consistently 
measured between BABAR and Belle, and only B(B° d -> D + D~ , D*±D*), C(B° d -> D + D~) and 
C(B d , B d — > D + D*~) have low consistency between BABAR and Belle. 

Our SM estimates predicted within the theoretical uncertainties of input parameters are 
given in the second columns of Table [3] and Table HI Theoretical predictions for the branch- 
ing ratios and the polarization fractions are given in Table El CPA predictions are given in 
Table HI All the branching ratios are above 10~ 4 order. The direct CPAs are expected to be 
quite small. All mixing- induced CPAs of B d decays are very large (about —0.7). There is an 
obvious signature of the mixing-induced CP violations in B° s — > D* + D~,D+Dl~ decays since 
S(B®,B® — > D* + D~) ^ —S(B®,B® — * D+D*~), which are consistent with the experimental 
measurements. In addition, for B d — > D* + D* , B~ — > D* D* and B® — > D* s + D* decays, 
the longitudinal and transverse polarization fractions can be precisely predicted, and are about 
~0.5 and ~0.1, respectively. Comparing present experimental data in Table [2] with the SM 
predictions in Table [3] and Table HI we can find that all measured quantities agree with the SM 
expectations within the error ranges except C(B d , B d — > D + D~) from Belle. 

We now turn to explore the RPV effects in B° d -> D M+ D ( * } ~ , B~ -> D ( * )0 D M ~ and B° s -> 
jj* + jj * (jgpayg fj^g mos t conservative existing experimental bounds are used in our analysis. 
We choose the averaged data, which have highly consistent measurements between BABAR and 
Belle (defined as a scale factor S < 1), and varied randomly within 2a ranges to constrain 
the RPV effects. The current experimental data and theoretical input parameters are not 
yet precise enough to set absolute bounds on the relative RPV couplings. We obtain the 
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Table 2: Experimental data for B® — > D l * )+ D {,) and B~ — > D* )0 D^* decays from BABAR and 
Belle. The branching ratios (£>) are in units of 10~ 4 . The scale factor S is defined in introduction 
part of Ref. 42|, and S > 1 often indicates that the measurements are inconsistent. 



Observable 


BABAR 


Belle 


Average 


S 


B(B°, -> 


D+D~) 


2.8 ±0.4 ±0.5 [9] 


1.97 ±0.20 ±0.20 [3J 


2.1 ±0.3 


1.2 


"-^V^d 




5.7±0.7±0.7 [9] 


11.7± 2.6±|;| [26] 


6.1 ± 1.5 


1.6 


"\ d 


D* + D*~) 


8.1 ±0.6 ±1.0 [9] 


8.1 ± 0.8 ± 1.1 [7] 


8.1 ±0.9 


< 1.0 




D°D~) 


3.8 ±0.6 ±0.5 [9] 


3.85 ±0.31 ±0.38 [10] 


3.8 ±0.4 


< 1.0 




D* a D~) 


6.3 ± 1.4 ± 1.0 [9] 








*^ V it 


D°D*-) 


3.6 ±0.5 ±0.4 [9] 


4.57 ±0.71 ±0.56 [11] 


3.9 ±0.5 


< 1.0 




D*°D*~) 


8.1 ±1.2 ±1.2 [9] 








" V d 


D+D-) 


-0.07 ±0.23 ±0.03 [4] 


-0.91 ±0.23 ±0.06 [3] 


-0.48 ±0.42 


2.5 


C(B°B° 


-f D* + D~) 


0.08 ±0.17 ±0.04 [4] 


-0.37 ±0.22 ±0.06 [5j 


-0.09 ±0.22 


1.6 


^ V^d i -"d 


-> D+D*-) 


0.00 ±0.17 ±0.03 [4] 


0.23 ±0.25 ±0.06 [5J 


0.07 ±0.14 


< i.o 


^ V d 


> D* + D*~) 


0.00 ±0.12 ±0.02 [4] 


-0.15 ±0.13 ±0.04 [8] 


-0.07 ±0.09 


< i.o 


"MJP v it 


-> D°D~) 


-0.13 ±0.14 ±0.02 [9] 


0.00 ±0.08 ±0.02 [10] 


-0.03 ±0.07 


< i.o 


•^CP \ u 


-> D*°D~) 


0.13 ±0.18 ±0.04 [9] 








•^CP \ u 


-> D°D*~) 


-0.06 ±0.13 ±0.02 [9] 


0.15 ±0.15 ±0.05 [11] 


0.03 ±0.10 


< i.o 


A+f n (B 


- -> D*°D*-) 


-0.15 ±0.11 ±0.02 [9J 










D+D~) 


-0.63 ±0.36 ±0.05 [4] 


-1.13±0.37±0.09 [3J 


-0.87 ±0.26 


< i.o 


S{B^B a 


->■ D*+D-) 


-0.62 ±0.21 ±0.03 [4] 


-0.55 ±0.39 ±0.12 [5J 


-0.61 ±0.19 


< i.o 


S{BlB[ 


-> D + D*~) 


-0.73 ±0.23 ±0.05 [4] 


-0.96 ±0.43 ±0.12 [5] 


-0.78 ±0.21 


< i.o 


S + (B° d ~ 


> D*+D*-) 


-0.76 ±0.16 ±0.04 [4] 


-0.96±0.25±g:ig [8] 


-0.81 ±0.14 


< i.o 


f±{B° d - 


■ D*+D*-) 


0.158 ±0.028 ±0.006 [4] 


0.125 ±0.043 ± 0.023 [8] 


0.150 ±0.025 


< i.o 




D*+D*-) 




0.57 ±0.08 ±0.02 [7J 







allowed scattering spaces of the RPV couplings A 2 3 2 A 212 and A^A^i as displayed in Fig. [TJ 
These survived parameter spaces are not in conflict with the above mentioned highly consistent 
experimental data in B% — > D ( * )+ D ( * > and B~ — > D* °D* decays. 

The squark exchange coupling A 2 3 2 A 212 contributes to all twelve B® — > D ( * )+ D* , B~ — ► 
D*°D* and B Q S — > D* + D * decay modes. The allowed space of A 2 3 2 A 2 \ 2 is shown in the left 
plot of Fig. [TJ Its magnitude | ^232^212 1 an d its RPV weak phase 0rpv have been constrained 
significantly. We obtain |A^ 2 A 2 ' 12 | G [0.14, 1.62] x 1CT 3 and RPV G [-75°, 84°]. The right plot of 
Fig. [TJ displays the allowed space of the RPV couplings A^ 23 A^ 21 due to slepton exchanges, which 
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Table 3: Theoretical predictions for CP-averaged branching ratios (in units of 10~ 4 ) and ratios of 
polarization (in units of 10 -2 ) in exclusive color-allowed b — ► ccd decays. The second column gives the 
SM predictions with the theoretical uncertainties of input parameters. The last two columns are the 
RPV MSSM predictions with different RPV couplings considering the input parameter uncertainties 
and experimental errors. 



Observable 


SM 


MSSMw/ A^ 2 A' 2 ' 12 


MSSMw/ Af 23 A' m 


B(B° d - 


D+D-) 


[2.35,4.15] 


[2.77,3.80] 


[2.77,4.39] 


B(3° d - 


D*+D-) 


[2.27,3.96] 


[2.87,4.22] 


[2.30,4.59] 


B(B° d - 


D + D*~) 


[2.56,5.04] 


[3.31,4.65] 




B(B° d - 


D* ± D T ) 


[4.84,8.95] 


[6.18,8.70] 


[5.21,8.84] 


B(B° d - 


D* + D*-) 


[6.21, 12.22] 


[7.05,8.90] 




B(B- - 


■ D°D~) 


[2.53,4.43] 


[3.00,4.04] 


[3.00,4.68] 


B(B~ - 


■ D*°D-) 


[2.42,4.27] 


[3.07,4.53] 


[2.25,4.75] 


B(B- - 


■ D°D*-) 


[2.73,5.42] 


[3.55,4.96] 




B(B~ - 


■ D*°D*-) 


[6.61,13.10] 


[7.54, 10.67] 




B(B° S -> 


DfD-) 


[2.33,4.20] 


[2.76,3.81] 


[2.77,4.48] 


B(B° S - 


D*+D-) 


[2.24,4.00] 


[2.81,4.25] 


[2.08,4.42] 


B(B° S - 


DfD*-) 


[2.54,5.03] 


[3.27,4.69] 




B(B° S -> 


D*+D*-) 


[6.15,12.10] 


[6.92, 10.04] 




urn - 


> D* + D*-) 


[52.40,52.97] 


[50.35,52.53] 




h(B- - 


■* D*°D*-) 


[52.43,53.02] 


[50.37,52.56] 




h(B° s - 


> D* + D*-) 


[52.56,53.16] 


[50.60,52.70] 




f±(B° d ~ 


* D* + D*-) 


[8.82,9.51] 


[10.00,12.94] 




h(B- - 


-> D*°D*-) 


[8.84,9.53] 


[10.03,12.97] 




h(B° s - 


* D*+D*~) 


[8.35,9.05] 


[9.50, 12.34] 





contributes only to six decay modes B® — > D { * >+ D , B~ —> D*° D and B® — > D { * )+ D . The 
magnitudes |A / i23 A^ 21 | have been limited within |A^ 23 A^ 21 | < 1.28 x 10~ 3 , and the corresponding 
RPV weak phase </>rpv for the range |A' i23 A^ 21 | < 0.4 x 10~ 3 is not constrained so much, however, 
the RPV weak phase for |A^ 23 A' i21 | G [0.4, 1.3] x 10~ 3 is very narrow. 

Using the constrained parameter spaces shown in Fig. one may predict the RPV effects 
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Table 4: Theoretical predictions for CPAs (in units of 10 2 ) in exclusive color-allowed b — ► ccd decays. 



Observable 


SM 


MSSM w/ A 2 '$ 2 A 2 ' 12 

/ _•>_. _ j_ _ 


MSSM w/ A'* 2 ,A' 21 




+ D + D-) 


[-78.00, -71.67] 


[-97.52,-52.66] 


[-99.83,-35.16] 


S(B°,,B°, - 


* D* + D-) 


[-70.40, -64.55] 


[-81.77,-32.17] 


[-98.01,-55.02] 


S(B°,,B , - 

K-^di d 


* D + D*-) 


[-72.17,-66.83] 


[-83.29,-36.15] 


[-98.30,-57.69] 


S + (B°,,B® 


D* + D*-) 


[-72.73, -67.77] 


[-95.18,-53.01] 




C(B°,,B°, - 

^\ di a 


> L>+L>-) 


[-6.03,-3.87] 


[-7.61,0.92] 


[-11.05,2.59] 




> D* + D-) 


[3.36, 13.83] 


[1.38,14.75] 


[-13.35,21.30] 


C(B°B°, - 


> D + D*-) 


[-14.44,-3.53] 


[-14.83,-1.45] 


[-21.00,13.28] 


C + (B°,,B°, 


£)*+£)*-) 


[-1.35,-1.03] 


[-1.83,0.20] 




A^UB- -s 


D°D-) 


[3.87,6.03] 


[-0.92,7.61] 


[-2.59,11.05] 


AnUB- -> 

■'MJP \ u 


D*°D~) 


[-1.15,-0.45] 


[-1.10,0.29] 


[-1.99,0.23] 


AnUB~ -S 
■'MJP \ u 


D°D*-) 


[1.03,1.35] 


[-0.40,1.42] 




App r (B~ 


-> D*°D*~) 


[1.03,1.35] 


[-0.20,1.83] 




A^{B° S - 


D+D-) 


[3.87,6.03] 


[-0.92,7.61] 


[-2.59,11.05] 


A&(B° S - 


D*+D-) 


[-1.15,-0.45] 


[-1.10,0.29] 


[-1.99,0.23] 


-> 


D+D*-) 


[1.03,1.35] 


[-0.40,1.42] 




A + c f v (B° - 


-> D*+D*-) 


[1.03,1.35] 


[-0.20,1.83] 






0.0 0.5 1.0 1.5 2.0 0.0 0.3 0.6 0.9 1.2 1.5 



l^ 2 3 2 r 2 1 2 K X10 " 3 ) l^ i2 3^ i2 lK X10 " 3 ) 

Figure 1: Allowed parameter spaces for relevant RPV couplings constrained by — > D W+ D* and 
B~ — > D ( ^° D ) , where </>rpv denotes the RPV weak phase. 
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on the other quantities which have not been measured yet or have less consistent measurements 
between BABAR and Belle. With the expressions for B, C, S, A^p, fh and f± at hand, 
we perform a scan on the input parameters and the constrained RPV couplings. Then we 
obtain the RPV MSSM predictions with different RPV coupling, whose numerical results are 
summarized in the last two columns of Table [3] and Table HI 

The contributions of A^A^ due to squark exchange are summarized in the third columns of 
Tableland Table HI In Tabled comparing with the SM predictions, we find A^A^ coupling 
could not affect all branching ratios much. Three Jl{B — > D*D*) and three f±(B — > D*D*) are 
slightly decreased and increased by A^A^ coupling, respectively, and their allowed ranges are 
scarcely magnified by this coupling. As given in Table HJ the A^A^ contributions could 
greatly enlarge the ranges of four S(B®,B® — » D ( * )+ D ( * ] ). The effects of Af^A^ cou- 
pling could extend a little bit the allowed regions of four C(B%,B° d -> £>(*)+£>(*)-) and eight 
Acp(B- -> D^°D^-,B° S -> dW+DW"), too. But this squark exchange coupling cannot 
explain the large C(B®, B d — > D + D ) from Belle. The predictions including slepton exchange 
couplings A^ 3 A^ 21 are listed in the last columns of Table [3] and Table HI The A^ 3 A^ 21 couplings 
do not give very big effects on the relevant branching ratios, but could significantly magnify 
the ranges of S(B®,B® — > D ( * )+ D ) from their SM predictions as well as extend the ranges 
of C(B° d ,B° -> £>(*)+£>-) and A^(B~ -> D^ D~,B S -> D^ + D~). The lower limits of 
C(B d , Bj — > D^ + D~) could be reduced by these slepton exchange couplings, too, but slepton 
exchange coupling effects are still not large enough to explain the large C(B d ,B d — » D + D ) 
from Belle. 

It is worth noting that our investigation of the color-allowed b —>■ ccd decays was motivated 
by the large direct CPA of B® — > D + D reported by Belle which has not been confirmed 
by BABAR and contradicted the SM prediction. Relative RPV couplings, constrained by all 
consistent measurements in B d — > D* + D w and B~ — > D*°D^* ) systems, could slightly 
enlarge the range of C(B d ,B d — > D + D ). Our RPV MSSM prediction for this observable is 
coincident with the BABAR measurement, but still cannot explain the Belle measurement. The 
unparticle interaction has positive effects on C(B d , B® — > D + D ) as obtained in Ref. [jjj, in 
which the author used only experimental constraints of B(B d — > D + D ). Note also that very 
large value of C(B d ,B d — > D + D ) could be obtained by unparticle interaction, however, with 
the sign opposite to the Belle measurement. 

For each RPV coupling product, we can present correlations of physical quantities within the 
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constrained parameter spaces displayed in Fig. [T] by the three-dimensional scatter plots. RPV 
coupling contributions to B d — > D { * )+ D M , B~ — > D* )0 D { * and — > D* )+ D { * decays are 
very similar to each other. So we will take an example for a few observables of B d — > D ( * )+ D *' 
decays to illustrate RPV coupling effects. Effects of RPV couplings A^A^ and A^ 23 A^ 21 on 
observables of B d — > D* + D* decays are shown in Fig. [2] and Fig. [HI respectively. 

In Fig. El we plot B, ft,, C and S as functions of A 2 3 2 A 212 . Three-dimensional scatter plot 
Fig. El (a) shows B(B d — > D + D~) correlated with |A 2 3 2 A 212 | and its phase 0rpv- We also give 
projections to three perpendicular planes, where the IA^A^I'^rpv plane displays the con- 
strained regions of A 2 3 2 A 212 , as the left plot of Fig. [TJ It is shown that B(B d — > D + D~) 
is little decreasing with | -^232-^212 1 011 the I ^232^212 1 plane. From the £>-0rpv plane, we 
see that B(B d — > D + D~) is not sensitive to 0rpv- All other branching ratios have simi- 
lar trends with A 2 3 2 A 212 coupling. From Fig. El (b-d), we can see fi(B d — > D* + D*~) and 
C(B° d ,B Q d — > D + D~,D* + D~) are not very sensitive with |A 2 3 2 A 212 | and 0rpv- RPV coupling 
A232A212 contributions to S(B d ,B d — > D + D~ , D* + D*~) are also very similar to each other. 
So we take an example for S(B d ,B d — > D + D~) shown in Fig. El (e) to illustrate the RPV 
coupling effects. S(B d ,B d — > D D ) is decreasing (increasing) with | -^232^212 1 when 0rpv > 
(0rpv < 0), and it is decreasing with 0rpv- S(B° d , B a d D* + D~, D + D*~) have totally different 
trends to S(B%, B° d -> D D~) with (A^A^I and 

0rpv, and we show only the squark exchange 

effects on S(B° d , B° d -> D* + D~) in Fig. El (f). 

Fig. [3] gives the effects of the slepton exchange couplings A^A^ in B d — > D ( >+ D decays. 
As displayed in Fig. [3] (a), B(B d — > D + D ) is not very sensitive with (A^A^il and has only 
small allowed values when |</>rpv| is small. Fig. [3] (b) shows that B(B d — > D* + D ) is decreas- 
ing with |A^23A'j 2 i| and is weakly sensitive to |0rpv|- Fig. [3J (c) exhibits the slepton exchange 
coupling effects on C(B d , B d — > D + D~), which is decreasing with |Aj£ 3 Aj 2 i| and has little sensi- 
tivity to 0rpv- Slepton exchange couplings have great effects on C(B d , B d — > D* + D~ , D + D*~) 
and S(B d ,B d — > D + D ,D + D ,D + D* ), and they have quite complex variational trends to 
|A-*2 3 A-2il and |0 RPV |- The effects of Af 23 A^ 21 couplings on C{B° d , B° d -> D*+D~) and S{B%, B° d -> 
D + D~ ,D* + D~) are shown in Fig. [3j(d-f). C{B G d ,B Q d -> D + D*~) has entirely different trends 
from C(5°, 5° -> D* + D~). S(B° d , B° d -»■ D + D*-) has a similar trends as 5(5°, 5° -> D* + D-). 
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3.2 Exclusive color-allowed b — ► ccs decays 



Exclusive color-allowed b — > ccs tree decays include B® — > D M+ D* , £>~ 



D* ° D { * } and 



D ( J )+ decay modes. Almost all branc hing ratios and one longitudinal polarization 



have been measured by Belle 49j , BAB AR 
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29], CLEO 
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collaborations. Their averaged values from Particle Data Group 42j] are listed as follows 
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33|, and ARGUS [34] 







= (7.4 ±0.7) x 10" 3 , 


B(B° d 


B(B° d - 


*D + D*-) 


= (7.6 ± 1.6) x 10" 3 , 


B(B~d 


B(B- 




= (10.3 ± 1.7) x 10~ 3 


, B(B 


B{B~ - 


-> D°D*~) 


= (7.8 ± 1.6) x 10" 3 , 


B(B~ 


B(B°s 


-> DtD;) 


= (ll±4)xl0~ 3 , 


B(B° S 


B(B° S - 




< 257 x 10^ 3 , 


h(B° d 



d* + d;) = (8.3 ± 1.1) x 10~ 3 , 

D* + D*~) = (17.9 ± 1.4) x 10~ 3 , 
-> D*°D^) = (8.4 ± 1.7) x 10" 3 , 
> D*°D* S ~) = (17.5 ± 2.3) x 10" 3 , 

DjD*-) < 121 x 10~ 3 , 

D* + D*-) = 0.52 ±0.05. 



(33) 



The SM predictions, in which the full theoretical uncertainties of input parameters are 
considered, are given in the second columns of Table [5] and Table El Theoretical predictions 
for the branching ratios and the polarization fractions are given in Table El Predicted CPAs 
are also given in Table El Compared with the experimental data, only the SM predictions of 
13(13® — > D* + D*~ , B~ — > D*°D*~) are slightly larger than the corresponding experimental data 
given in Eq. (133|) . and all the other branching ratios are consistent with the data within la error 
level. For the color-allowed b — > ccs decays the penguin effects are doubly Cabibbo-suppressed 
and, therefore, play a significantly less pronounced role in CPAs. These CPAs have not been 
measured yet. We obtain that all CPAs are expected to be very small (about 10~ 3 or 10 -4 
order) in the SM except S(B° S ,B° S -> D* + D*~) and C(B° S ,B° S -> D* + Dj , D+ D*~) . There is 
no obvious signature of CP violation in B s — > D* decays since C(B®,B® — > D* + D~) m 
-C(B°, 4° - DfDr ) and S(B° S , B° s -> D*+Dj) » -S(B° S , B° s -> D+D* s ~). 

There are two RPV coupling products, A^A^i and A^A^, contributing to these exclusive 
b — > ccs decay modes at tree level. We use the experimental data listed in Eq. (1331) to constrain 
the RPV coupling products, and the allowed spaces are shown in Fig. HI The coupling A^A^ 
due to squark exchange contributes to all twelve relative decay modes. The allowed space of 
A231A221 is shown in the left plot of Fig. HI The slepton exchange couplings X'^ 3 \' i2 2 contribute 
to six B® — > D ( * )+ D s , B~ — > D { * )0 D s and B® — > D ( * )+ D s decays, and the constrained space is 
displayed in the right plot of Fig. HI From Fig. HI we find both moduli of RPV couplings have 
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Table 5: Theoretical predictions for CP averaged B (in units of 10 4 ) and polarization fractions (in 
units of 10 -2 ) of exclusive color-allowed b — > ccs decays in the SM and the RPV MSSM. 



Observable 




MoblVl W/ ^23x^221 


JVlbbJVl w/ A i2 3A i2 2 


\ o. 


* D + D s) 


rft 7n m eel 

|o. / u, iu.ooj 


[o.oo, /.oyj 


rc /in o onl 
[0.4Z, o.oUJ 


B(B° d - 


> D* + D;) 


[0. / U, 1U.40J 


\a A7 n /i nl 

[o.4 < , y.4yj 


rc i c n Qnl 
[o.lo, y.oUJ 


B(B° d - 


> D+D*-) 


TV QO 1 Q OOl 

[i.oz, lo.ZZJ 


re on 1 n onl 

[o.yu, lu.zyj 




B(B° d - 


* D* + D*~) 


[in 17 Q/l /I Ol 

[ly.Z / , o4.4ZJ 


M s co on vnl 
[lo.oy, ZU. / UJ 




B(B~ - 




rv 01 1 1 /i qi 
[1.21, 11.4oJ 


[fi nn c. 1 ol 
[o.yu, o. 1ZJ 


\r on o /iol 

[o.yu, y.4yj 


B(B- - 


-> D*°D~) 

S J 


[/ .1 1 , 11.Z4J 


re n c; 1 n 1 71 
[o.yo, 1U.1 i\ 


\R RA O OKI 

[o.o4, y.yoj 


B(B~ - 


-> D°D*~) 


on i /i OTl 
[f.oy, 14./ (\ 


\i a o 11 nnl 
[/.4o, il.UUJ 




B(B- - 


-> D*°D*-) 


Ton q7 n^l 
[21). I , o / .UoJ 


Ti n nn oo 1 nl 

[iy.yy, zz.iuj 




B(B° S - 




[o.oo, 1U. / Zj 


q£ v vi i 

[O.OO, 1 . 1 1J 


oo o ncl 

[O.zo, y.ubj 


B(B° S - 




\R A R ~[ C\ A A~\ 

[0.4o, 1U.44J 


1 n /i £l 

[o.oi, y.4oj 


\r no n o^l 

[o.uz, y.zoj 


B(B° S - 




Tv no 1 o nvl 
[ f .Uo, Iz.y /J 


Tv nn 1 n a nl 
[ / .UU, 1U.4UJ 




B(B° S - 


> D* s + D*~) 

S S J 


[lo.o4, oo.ooj 


\~\ q aq on nol 

[lo.4o, zu.yoj 




h{B° d ~ 


D* + D*~) 


rcn op; cn n 1 1 

[oU.Zo, ou.yij 


\ AQ AC CI 1 Ql 

[4o.4o, ol.loj 




fr(B~ 




rcn oq cn n/il 
[oU.Zo, 0(J.y4J 


\ AQ Ad C1 1 ^1 

[4o.4y, ol.loj 




h(B° s - 




[50.40,51.10] 


[48.71, 51.30] 




f±(B° d - 


D* + D*~) 


[8.85,9.55] 


[8.15,12.85] 




f±(B~ 


D*°D*-) 


[8.87,9.57] 


[8.17,12.88] 




f±(B° s - 




[8.38,9.07] 


[7.71,12.23] 





been limited as lA^A^il < 8.05 x 10~ 3 and |A-2 3 A / i22 | < 5.05 x 10 -3 . Their RPV weak phases 
are not constrained much when their magnitudes are less than about 1 x 10~ 3 . 

Next, using the constrained parameter spaces shown in Fig. HJ we are going to predict 
RPV effects on the observables which have not been measured yet. We summarize RPV MSSM 
predictions with two separate RPV coupling contributions in the last two columns of Table 
and Table El 

The contributions of A^A^i coupling due to squark exchange are summarized in the third 
columns of Table E] and Table El In Table El we find that the ranges of all branching ratios 
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Table 6: Theoretical predictions for CPAs (in units of 10 2 ) of exclusive color-allowed b — ► ccs decays 
in the SM and the RPV MSSM. 



Observable 


SM 


MSSM w7 X'^X'L, 


MSSM w7 A'noA'no 


J&{B° d - 


D+D-) 


[_n 34 _n 221 


[— 3 06 2 581 


[—8 42 7 941 




D*+D-) 


[n 03 061 


[— 32 361 


[_0 98 1 071 


A&(B° d -> 


D+Dt) 


[_n 07 _n 061 


[— 51 441 




A + c ^(B° d - 


-> D* + D*-) 


[_o 07 —0 061 


[_0 69 561 




A d SUB- - 


• D°D~) 


r_o 34 —0 221 


[— 3 06 2 581 


[—8 42 7 941 


A d SUB- - 


D*°D;) 


[0 03 061 


1— 32 361 


[_0 98 1 071 


Acp(B~ — 


D°D*-) 


r_o 07 —0 061 


[_0 510 581 




A + c f v (B- - 
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Figure 4: Allowed parameter spaces for relevant RPV coupling products constrained by the measure- 
ments of exclusive color-allowed b — > ccs decays listed in Eq. (|33|). 
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are shrunk by A^A^i coupling and the experimental constraints. Especially, A^A^i coupling 
effects could reduce the range of B(B® — > D* s + D* s ). However, the allowed ranges of three 
/l(-S(s) — >• D* S )D*) and three /±(-B( s ) — > D* S ^D* S ) are enlarged by A^A^ coupling. In Table El 
we can see A^A^i coupling does not affect C(B Q S — > D* s + D s , D* s + D s ) much. 

Meanwhile, RPV coupling effects could remarkably enlarge the allowed ranges of the other 
direct CPAs (about 10 times). Unfortunately, they are still too small to be measured at 
presently available experiments. It is interesting to note that mixing-induced CPAs of B s 
decays are greatly affected by A^A^i coupling. For an example, \S{B° S — > D* + D* )| could 
be increased to ~ 50% and their signs could be changed by the squark exchange coupling. 

The contributions of A'^A^ due to the slepton exchanges are listed in the last columns 
of Table El and Table El From Table El we find the ranges of all branching ratios are shrunk 
by A'j 2 3A^22 coupling and the experimental constraints. The last columns of Table El show 
that A^A^ coupling could enlarge the ranges of all the CPAs. Particularly, A^A^ coupling 
could change the predicted S{B° S ,B® — > D* D s ) significantly from quite narrow SM ranges to 
[—0.6, 0.6] or [—1, 1]. The upper limits of \C(B®, B° s — > D* s + D s , D^D* a )| are increased a lot by 
A / i2 3 A / i22 couplings. 

Since RPV contributions to physical observables are also very similar in B® — > D w+ D* , 
B~ — > D*°D* and B Q S — > D* + D* systems, we show only a few observables of B s decays as 
examples. Fig. Eland Fig. El show the variational trends in some observables with the A^A^i 
and A / i 2 3 A / i22 couplings, respectively. 

First, we will elucidate the information implied in Fig. El From Fig. El (a), we find 
B(B° S — > D* + D*~) is not changed much by (A^A^il, and could have only small value when 
0rpv is not too large. As shown in Fig. El(b), f±(B® — > D* + D*~) is increasing with (A^A^I 
and also could have small value when 0rpv is small. From Fig. El (c-d), we find \S(B®,B® — > 
D + S D S ,D* S + D S )| are all rapidly increasing with | -^231^221 1 an d could be very large at the large 
values of (A^A^il and |0rpv|, furthermore, the RPV weak phase has opposite effects between 
\S(B° S ,B° S - D + 8 D~)\ and \S{B° S ,B° S -> D\ + D~)\. A^A^ coupling effects on S{B° S ,B° S -> 
D+D*~) (S(B° S ,B° S -> D* s + D*~)) are similar as ones on S(B° S ,B° S -> D* + D~) (S(B° S ,B° S -> 
D+Dj)). So any measurement of S(B®, B° s — ► D { * )+ D ( * ] ) in the future will strongly constrain 
the magnitude and RPV weak phase of A^A^i coupling, and then other mixing- induced CPAs 
will be more accurately predicted as indicated by Fig. El (c-d). As shown in Fig. El (e), 
C(B° S ,B° S -> D + S D~) has similar trends as S(B° S ,B° S -> D + a D~) with (A^A^I and RPV , 
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Figure 5: Effects of RPV coupling A^A^ in 5° -> D ( * )+ D { * ] decays, where lA'^A'^l is in units 
of 10 -3 , B in units of 10 -4 , and Jt denotes the transverse polarization fraction f±. 




Figure 6: Effects of A / i 2 3 A^ 22 in B® — > L> s * + ZA, decays, where lA^A^I are in units of 10 3 , and B in 
units of 10~ 4 . 
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however, A 231 A 221 coupling effects on the former are much smaller than the effects on the latter. 
Acp(B U7 d — > DD S ,D*D S , D* D* s ) and — Aq^{B u ^ — > DD* S ) have the same variational trends 
with X'^X'^i as C(B° S ,B° S -> D + ,D~) has. C(B° S ,B° S -> D* s + D~ , D* s + D~ ) are not affected much 
by A^Ag-j! coupling, and we show C(B°, £° -> £>*£>;) in Fi §- El (f) as an example. 

Fig. [6] illustrates A / i 2 3 A / i22 contributions to the CPAs of B® — > D* + D s . As displayed in Fig. 
El (a-b), S(B° S , B° s - are very sensitive to Af 23 A' i22 couplings. \S(B° S , B° s -> 

are strongly increasing with |A^ 23 A' i22 |, and they could reach extremum at |0rpv| ~ 120°. The 
A- 23 A- 22 coupling effects on S(B®, B° s — > D + S D*~) are same as the ones on S(B®, B° s — > D* s + D s ). 
Fig. M (c) shows that C(B®,B® — > D^D S ) are also very sensitive to |A' i23 A^ 22 | and 0rpv, 
but it is still too small to be measured in near future. In addition, A^>(B U) d — > D { * ] D s ) 
are affected much by A^ 23 A^ 22 couplings, and just RPV MSSM predictions of these quantities 
are very small. A^ 23 A^ 22 couplings could have similar impacts on C(Bg,B® — > D s + D s ) and 
-C(B° S ,B° S -> D* s + D~). We give these coupling effects on C{B° S ,B° S -> £>* +£) 7) in Fi S- El (d), 
which shows C(B° S ,B° S -> D* s D~) is increasing with |A' i23 A / i22 |, and could have large value at 
large |0 RPV |- 

4 Summary 

We have studied the twenty-four double charm decays i?^ — > D* + D* , 5,7 — > D*°D* and 
5° -> Dl* 1 D ( * } ~ in the RPV MSSM. We have treated these decays in the naive factorization 

s s (s) ^ 

and removed the known k 2 ambiguities in the penguin contributions via b — > qg*{^*) — > qcc 
by calculating its hard kernel b —* c + D q . Considering the theoretical uncertainties and the 
experimental error-bars, we have obtained fairly constrained parameter spaces of RPV couplings 
from the present experimental data, which have quite highly consistent measurements among 
the relative collaborations. Furthermore, using the constrained RPV coupling parameter spaces, 
we have predicted the RPV effects on the branching ratios, the CPAs and the polarization 
fractions, which have not been measured or have not been well measured yet. 

The investigation of exclusive color-allowed b — > ccd decays is motivated by the large direct 
CPA C(Bj, B® — > D + D~) reported by Belle, which has not been confirmed by BABAR yet and 
contradicted the SM prediction. Using the most conservative experimental bounds from B® — > 
D* + D* and B~ — > D { * >0 D* systems (choose only twelve highly consistent measurements 
between BABAR and Belle), we have first obtained quite strong constraints on the involved 
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RPV couplings A^A^ and A^ 23 A^ 21 from b — > ccd transition, due to squark exchange and 
slepton exchanges, respectively. Then, using the constrained RPV coupling parameter spaces, 
we have predicted the RPV effects on C(B d ,B d — > D + D~) and other observables, which have 
less consistent measurements or have not been measured yet. We have found that the lower limit 
of C(B° d ,B° d -> D^ + D~) could be slightly reduced by the RPV couplings. Our RPV MSSM 
prediction of C(B d ,B d — > D + D ) is consistent with BABAR measurement within la error 
level, but cannot explain the corresponding Belle experimental data within 3a level. We have 
also found that the contributions of A^ 2 ^2i2 an( ^ ^23^21 cann ot affect the relevant branching 
ratios much. A^A^ or AJ^A^i contributions could greatly enlarge the ranges of the relevant 
mixing- induced CPAs S(B d ,B d — > D { * )+ D M ) from their SM predictions, and these quantities 
are very sensitive to the moduli and weak phases of A^A^ and A^A^i- So more accurate 
measurements of S(B d ,B d — > D ( ' )+ D ( * > ) in the future will much more strongly constrain 
these RPV couplings, and then mixing-induced CPAs can be more accurately predicted as 
well. Effects of A^A^ coupling could slightly extend the allowed regions of four C(B d , B d — > 
£)(*)+£)(*)-) an d eight A^ r P (B- -> L>« £>(*)-, 5° -> DW+DW-). C(B° d ,B^ -> D~) are 
also sensitive to the slepton exchange couplings A^ 3 A^ 2 i, and their signs could be changed 
by these couplings. Additionally, three /l(-B( s ) — » D* S )D*) and three f±(B^ — > Df^D*) 
are decreased and increased by A^A^ coupling, respectively, and their allowed ranges are 
magnified by these couplings. 

For B° d -> D w+ ££ )_ , B~ -> D { * )0 D { * } ~ and 5° -> D? + D™~ decays, RPV couplings 
A231A221 an d A^23A^22 contribute to these decay modes. We have found A'^A^i coupling effects 
could apparently suppress the upper limit of B(B® —>■ D* s + D* s ), and could slightly enlarge the 
allowed ranges of three /l(-B( s ) — > D*^D*) and three f±(B^ — > Df^D*), nevertheless these 
quantities are not very sensitive to the changes of | A^A^ | and 0rpv- C(B^ — > D* s + D s , D^D* S ) 
are not evidently affected by the squark exchange A^A^i coupling, and their upper limits are 
increased a lot by the slepton exchange A'^A'^ couplings. RPV couplings A^A^i and A^A^ 
could greatly enlarge all other CP asymmetries, which are also very sensitive to the relevant 
RPV couplings. However, the direct CPAs are still too small to be measured soon. We could 
explore RPV MSSM effects from the mixing-induced CPAs of B s decays. 

With the large amount of B decay data from BABAR and Belle, especially from LHCb in 
the near future, measurements of previously known observables will become more precise and 
many unobserved observables will be also measured. From the comparison of our predictions 
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in Figs. [2H3 and Figs [MS with near future experiments, one will obtain more stringent bounds 
on the product combinations of the RPV couplings. On the other hand, the RPV MSSM 
predictions of other decays will become more precise by the more stringent bounds on the RPV 
couplings. The results in this paper could be useful for probing RPV MSSM effects, and will 
correlate with searches for direct supersymmetry signals at future experiments, for example, 
the LHC. 
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